Thermalization at lowest energies? 
A view from a transport model 
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Abstract 

Using the Isospin Quantum Molecular Dynamics (IQMD) model we 
analyzed the production of pions and kaons in the energy range of 1-2 
AGeV in order to study the question why thermal models could achieve 
a successful description. For this purpose we study the variation of pion 
and kaon yields using different elementary cross sections. We show that 
several ratios appear to be rather robust versus their variations. 

1 Introduction 

In the last decades heavy ion experiments have studied various combinations 
of heavy ion collisions. Besides the understanding of dynamical observables 
the production of particles have been of a special focus. Recent experiments 
at the SIS accelerator at GSI have us provided with high-quality data of the 
production of pions [T] and of kaons [51 [3] . The yields of particle ratios have 
been shown of large interest especially for the application of thermal models [4] . 
Several ratios of measured yields could be interpreted by the use of law-of-mass 
relations [5J |5] . That idea was later on supported by a detailed analysis from 
transport models [5] . 

We use the Isospin Quantum Molecular Dynamics model (IQMD) [51 [TU] in 
order to study the sensitivity of those ratios on the variation of several cross 
sections. If the system is near a thermal equilibrium, the influence of the cross 
sections is negligible. The dominant cross section in the energy range under 
study is besides the elastic nucleon-nucleon cross section the inelastic reaction 
of two nucleons forming a A and the corresponding inverse reaction NN <-> NA. 
The produced A may also decay into a nucleon-pion pair and the produced pion 
can again be absorbed by a nucleon forming a A: A <H> Ntt. 

Figure [JJ presents on the left side the parameterizations of the total pp 
(dashed) and pn cross sections (dotted line) and their respective inelastic cross 
sections forming a A (pp full line, pn thick dotted line). The differences between 
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Figure 1: Energy dependence of different cross sections used in IQMD. 

the total and the inelastic cross section is given by the elastic cross section. It 
should be noted that the cross sections of the absorption of a A are related to 
their production cross sections via detailed balance. The right hand side shows 
the cross sections for the production of a A by absorption of a pion. The in- 
verse reaction, the production of a pion by decay of a A is governed by its decay 
width Ta (around 120 MeV in the pole). In the following we use these parame- 
terizations fitted to experimental data of free nucleon-nucleon and pion-nucleon 
reactions and scale them by a global factor. This has been done either by scal- 
ing all cross sections or only by scaling a dedicated cross section. It should be 
noted that the scaling of the cross section for NN — > NA automatically scales 
the inverse reaction with the same factor in order to assure detailed balance. 
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Figure 2: Effect of the scaling of the cross sections on the isotropy ratio of the 
momentum space distribution in central Au+Au collisions. 
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2 Isotropy and pion production 



The first question is whether a dynamical equilibration in nuclear phase space 
has been reached or not. One criterion is an isotropy of the momentum dis- 
tribution in the center of mass. We present in Fig. [5] the ratio of transverse 
energy divided by twice the longitudinal energy, E t /(2Ei), as a function a mul- 
tiplicative scaling factor to the nucleon cross sections. We have chosen central 
collisions of a heavy system (Au+Au, b=Q fm at 1.5 AGeV) which has a high 
degree of stopping. Smaller systems do reach only lower values of this ratio. 
With the standard cross section we do not reach unity (green thick line) , and a 
much higher cross section is needed to approach this value. We find a significant 
influence of the scaling of the inelastic NN cross section (dashed line) which is 
related to the loss of kinetic energy when producing a A. An important part 
is due to the elastic cross section, while the pion-induced cross section (dashed 
line) does not show any sensitivity. No saturation is seen when further increas- 
ing the cross sections which would also be a possible signature of an equilibrium. 
We can thus conclude that a global dynamical equilibrium is not reached. 
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Figure 3: Left:The total pion yield in central Au+Au collisions at 1.5 AGeV as 
a function of the scaling parameter to the cross sections. Right: The total pion 
yield as function of the system size A in central A+A collisions at 1.5 AGeV. 
Different scaling parameters to the inelastic cross sections have been employed 

Next we study the influence of the cross sections on the pion yield, depicted 
on the left side of Fig. [3l and no saturation can be seen: The increase of the 
inelastic NN cross sections (NN — > NA and its inverse reaction NA —> NN, 
dashed line) does not enhance the pion yield significantly. Here we seem to reach 
a balancing of both channels. On the other hand the increase of the Nir — > A 
cross section (dotted line, its reverse reaction is governed by the A decay width 
Ta) decreases the pion yield. This reaction reinserts the A and enhances the 
probability of TV A — >• NN, while the decay A — >• Nir diminishes it. As a results, 
we change the balance point of the NN <->■ NA system. When applying the 
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scaling factor to all cross sections (full line), we obtain a maximum of the pion 
yield for a scaling factor of 1, i.e. for the present standard parametrization. 

Let us now study the influence of this scaling parameter to the system-size 
dependence as shown on the right side of Fig. [3J The pion yield has been 
divided by the participant number, which is equivalent to the system size for 
6 = 0. Only the standard parametrization (factor 1, full line) assures a scaling of 
the pion number with the system size as confirmed by experiment (see e.g. [1]), 
while other scaling factors results in increasing or decreasing curves. This special 
situation comforts the application of thermal models at these low energies, since 
the pion yield can now be described as an extensive variable. 





Figure 4: Comparison of the excitation function of the pion yield in central 
Au+Au collisions with FOPI data (top left) and of the K + (top right) and K~ 
(bottom) yield in with experimental inclusive data of the KaoS collaboration. 
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3 Production of K + and K 



Before continuing we would like to state that IQMD calculations are able to 
explain the experimental pion and kaon excitation functions within about 10- 
20% precision, as it is demonstrated in Fig. [4] The experimental data stems 
from the FOPI Collaboration for the pion [1] and for the K + and K~ from 
KaoS Collaboration [3]. 

The kaon production at this energy range is dominated by two-step processes 
like NiN 2 -> N 3 A -> NYK [TTJ [H]. This chain of reactions is sensitive 

to the yield of As and thus to a(NN •*-» NA). Since the production occurs in 
two steps, the production probability depends also on the density reached in the 
system. Therefore, the kaon yield is related to the degree of stopping indicated 
in Fig. [21 Furthermore, the kaon yield depends on the life time of a A, i.e. the 
duration when it is available for a collision. 
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Figure 5: Influence of a scaling factor to the A-decay width Ta on the final pion 
and kaon yield (left) and on the contributions to the kaon yield (right) 

The influence of the life time of the A on the kaon yield is demonstrated 
in Fig. [5j left, applying an multiplicative factor the decay width Ta which is 
about 120 MeV at the pole. A factor less than 1 decreases the width and 
enhances the life time. A larger factor reduces the life time and thus penalizes 
the reaction A A — > NN. the influence of the life time is partially balanced 
by reabsorption of the produced pion via the Nn — > A reaction as can be seen 
on the left side of Fig. [5] where we show the final yield of pions (dotted line, 
divide by a factor 400 to compare with the kaons). Only for drastic changes 
of the decay width the yield changes significantly. The yield of the kaons (full 
line) changes quite moderately. The reason can be seen on the right side of 
Fig. El The dominant channel AiA 2 A A and A 3 A -> NYK (dashed line) 
is sensitive to the A lifetime and drops for larger values of Ta ■ However, this 
effect is to a large extend balanced by the appearance of another production 
chain AiA 2 ->• AA, A -> Att, A 3 tt -> YK (dotted line). Only a small 
sensitivity of the total yield (full line) remains. 
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Figure 6: Excitation function of the K/ir ratio in Au+Au collisions for different 
scaling factors applied to NN <H> NA and irN — > A at the same time. 



The fact that the kaon production depends strongly on a scaling of the 
cross sections while the pions are only effected moderately explains the strong 
influence of a scaling to the ratio of K + /ir shown in Fig. |S1 The rising trend of 
the curves is similar, only the absolut values vary with the scaling factor. We 
would like to stress that the cross sections for the dominant chain including a 
NA — > NYK are experimentally unknown. Therefore, the absolute yield of the 
kaons contains a large uncertainty. Therefore, any conclusion from the absolute 
values of the K + /tt ratios in this energy range is quite difficult to establish. 
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Figure 7: Excitation function of the K~ /K + ratio in Au+Au collisions for 
different scaling factors to the cross sections and comparison of the K~/K + 
ratio with the tt /A ratio. 

The production of antikaons is dominated by charge exchange of the hy- 
perons: irY <H- NK. Since the yield of the hyperons is related to the yield of 



6 



kaons via associate production BB — > NYK (which has a much lower threshold 
than the direct production BB — »■ NNKK), there exist a direct link between 
the yield of kaons and that of antikaons via the hyperons. Here we find first 
indications for the validity of a law of mass action, as it has been reported in 
[51 [3] ■ This idea was supported by IQMD calculations [7J which showed, that 
the yield of antikaons is rather insensitive on scaling factors to the cross section 
of the reaction irY o NK if we apply it to both directions. This observation 
is interpreted as a dynamical balance between pions and hyperons on one side 
and nucleons and antikaons on the other side. 

A scaling of other cross sections (like iVA or Nir) does not effect that link 
between kaons and antikaons. This is demonstrated on the left side of Fig. 
where a scaling of the cross sections only causes a small effect on the ratio of 
K~ / K + . The other parameter governing this ratio is the pion yield, as shown 
on the right side of Fig. [7J experimental data of the ratio K~ /K + show the 
same excitation function as the experimental pion yield, scaled with a factor of 
the participant size [5]. 




Figure 8: Excitation function of the (K /K + )/(ir /tt + ) double ratio in Au+Au 
collisions for different scaling factors to the cross sections and in comparison to 
THERMUS results [13] 

The values of this ratio correspond to the values obtained by thermal models 
like THERMUS Qj]. Figure[8]shows the comparison of the ratio (K~ /K+)/(n- /tt+) 
between THERMUS (triangles, compiled from K~ /tt~ and K + /ir + taken from 
[13] ) with IQMD calculations (lines). Since the ratio 7r~/7r + is mostly sensitive 
to the isospin composition of projectile and target, the excitation function of 
the ratio (K~ / K + )/ (ir~ /ir + ) corresponds mainly to the ratio K~/K + . This 
ratio turn out to be independent from a scaling of the cross sections and the 
values are comparable to the values obtained with THERMUS. This observable 
demonstrates well the validity of the law of mass action. Due to this relation, a 
description with thermal models appears applicable. 
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4 Conclusion 



We have used the IQMD model to study the influence of a scaling of cross 
sections on the yield of pions, kaons and antikaons. The pion yield is slightly 
effected by a scaling of the cross sections and yields a maximum for the standard 
paramctrization obtained from the free experimental cross sections. The scaling 
of the pion yield with the participant number can only be assured when no 
scaling of the cross sections is done. The ratio of K/w depends strongly on the 
scaling of the cross sections but its absolute value suffer from the uncertainty 
in unknown production channels of the kaon. The ratio K~ / K + reflects nicely 
the properties of the law of mass action. Consequently, thermal models are 
applicable for some but not all ratios is this energy range. 
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